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Abstract 



The poles of the quantum scattering matrix (S'-matrix) in the complex mo- 
mentum plane have been studied extensively. Bound states give rise to S-matrix 
f3 poles, and other poles correspond to non-normalizable anti-bound, resonance and 

^ anti-resonance states. They describe important physics, but their locations can 

qh be difficult to find. In pioneering work, Nussenzveig performed the analysis for a 

~~' square well/ wall, and plotted the flow of the poles as the potential depth/height 

_^ varied. More than fifty years later, however, little has been done in the way 

of direct generalization of those results. We point out that today we can find 
such poles easily and efficiently, using numerical techniques and widely available 
software. We study the poles of the scattering matrix for the simplest piece- 
wise flat potentials, with one and two adjacent (non-zero) pieces. For the flnite 
O well/wall the flow of the poles as a function of the depth/height recovers the 

results of Nussenzveig. We then analyze the flow for a potential with two inde- 
pendent parts that can be attractive or repulsive, the two-piece potential. These 
^ examples provide some insight into the complicated behavior of the resonance, 

anti-resonance and anti-bound poles. 
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1 Introduction 



The scattering matrix has a rich, interesting history. It has even, in the past, been 
postulated to provide a fundamental physical viewpoint [1]. During the first half of 
the 20th century quantum field theory was 'plagued' by infinities. Many were skep- 
tical of the prospects of quantum field theory to explain physical reality. As a result 
the ^-matrix became of central importance and was extensively studied. Once renor- 
malization techniques resolved the problems with the infinities and more and more 
elementary particle phenomena were successfully explained using quantum field the- 
ory, the ^-matrix lost its fundamental role. It remains important in high energy physics 
and other fields, however, and is still an interesting object of study. For one, it pro- 
vides a useful perspective on the quantum physics of local interactions. At a more 
fundamental level, the poles of the S-matrix provide a unified description of stable 
and decaying states. Also, they are actively studied by mathematicians in relation to 
various aspects of spectral theory |2]. 

It is well known that bound states correspond to poles of the scattering matrix. 
However, the ^-matrix has other poles that are not associated with normalizable states, 
but still encode important and interesting physics P,!?]. In this paper we consider the 
positions of all the poles of the ^-matrix for certain elementary potentials, and their 
flow under deformation of the potential parameters. 

What is the physical significance of the S'-matrix poles that are not associated 
with bound states? In the late 1920's Gamow proposed an explanation of a-decay 
in terms of solutions of Schrodinger's stationary equation with complex eigenvalues, 
that satisfied a purely outgoing boundary condition, i.e. far enough from the origin 
the solutions were outgoing plane wave. These solutions could be thought of as wave 
functions, i.e. as the position representations of certain generalized state vectors. Those 
same Gamov vectors turned out to correspond to the poles of the S'-matrix and the 
residues of the propagator. What Gamow sought to apply to a-decay was in general 
a way of describing bound and quasi-stable states that emphasizes their similarities. 
The concept of a quasi-stable state is a fundamental one, and so it has been applied 
extensively, and in all areas of physics 

The physical effects of the S'-matrix poles unrelated to bound states are undeniable 
in scattering (see [9j, e.g.). Strictly speaking, however, the quasi-stable states associ- 
ated with those poles, the resonance and anti-bound states, are not true statesj^ The 
non-normalizability of their wave functions is one marked difference from the physical 
bound states. This non-unifying characteristic, however, can be tamed somewhat, in 
a mathematical way, by regularizing the integrations]^ interpreting the probability in 
a time-dependent setting pjj, and/or continuing to complex potentials (see [12j, e.g.). 
Upon continuation to complex potentials, one can also relate the different poles to each 
other (bound state poles to anti-bound state poles, e.g.) |13] . 

We should note that, in spite of their non-normalizability, the wave functions asso- 

""^Of course, that continues to this day. To mention some recent applications, they have been used 
to calculate tunneling ionization rates |S], to understand the phenomenon of diffraction in time 
to describe the decay of cold atoms in (quasi-) one-dimensional traps |3, and are directly relevant to 
recent condensed-matter experiments [S]. 

^For brevity and simplicity, however, we will continue to refer to them as states, rather than as 
'states' or virtual states, and rely on the context to make the distinction. See PU] for a discussion. 

^This was first done by Zel'dovich; see 
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dated with resonance and anti-bound states can be useful. After Gamow's work, they 
figured in expansions of the S'-matrix ([2], e.g.) and propagators ([ISj, e.g.) and other 
physical quantities. As a result, various ways of determining and approximating them 
have been developed (see [16j, e.g.). We will not consider them here, however. 

We focus on the S'-matrix. Remarkably, it can be determined by its poles-their 
locations and residues |3]. The pole locations, and their flow on deformation of the 
potential, are therefore clearly of interest. In the simplest case of a square well or wall 
(barrier), the movements of the poles in the complex momentum plane were studied in 
[T?j|^ The depth of a simple square well and the height of a square barrier were varied to 
change the position of the S'-matrix poles and generate their flow. More recently, certain 
of Nussenzveig's results were reproduced using simple, elegant graphical methods in 
[18| . and different poles were related by complex continuation in |13) . 

More than fifty years after Nussenzveig's work, however, little has been done in the 
way of direct generalization. Even for the next-simplest potentials, numerical methods 
must be used to find the S-matrix poles. But today, it is not difficult to perform 
the analysis in [T7] on a personal computer, using widely-available software (such as 
Mathematica, Maple, MATLAB, e.g.). To illustrate this point, we will reproduce those 
results here and generalize them to the next-simplest cases. We will plot the trajectories 
of all poles for two types of elementary piece- wise fiat potentials-attractive/repulsive 
(following |17| ) and a two-piece combination of both. For short, we will refer to the 
potentials as one-piece, or wall or well or wall/well, and two-piece, or well+wall, etc. 

The above potentials are interesting because they are simple, but still generic. Their 
simplicity permits a relatively easy numerical treatment of the pole structure and a 
complete description of the pole positions in the complex momentum plane, with a flow 
that is a function of only a few parameters. Most importantly, unlike certain special 
potentials such as the Dirac delta and the Coulomb potential, they are generic enough 
to exhibit all the possible families of poles. 

We should mention that another way to generalize the Nussenzveig flow is to study 
potentials for which the S'-matrix can be found exactly (see [19j, e.g.). Such solvable 
potentials are special, however, and so may show non-generic properties. Here, we 
prefer to generalize the wall/well in a very direct, simple way that we hope will not 
introduce any special features into the flow. 



Let a spin- zero, massive particle move on the real line with coordinate x. The Hamilto- 
nian of the system H = /2m + V{x) has a potential V{x) that vanishes fast enough]^ 
at infinity for scattering to be possible. The stationary Schrodinger equation 



determines the wave functions of the stationary states and their energies. Let us define 
the momentum variable (wave number) k so that k"^ = 2mE/h'^, as it will have an 

*In fact Nussenzveig did more-a spherically symmetric rectangular potential well/wall for zero 
angular momentum was worked out. Some results were also given for higher angular momenta. 

^We will adopt the same convention as [1], ch. 3, i.e. faster than 1/x. That guarantees asymptotic 
solutions of the type e='=*'^^. 



2 Background 
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Figure 1: Potentials studied. The dashed line shows the square well/wall potential with 
variable depth/height U > 0/U < 0. The solid line represents a 2-piece flat potential with an 
attractive part and repulsive part, the well+wall. 
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Figure 2: Types of S -matrix poles in a typical configuration. Different types of poles 
correspond to wave functions with different boundary conditions. 
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important role in this work. We will set m = h = 1 in anticipation of the numerical 
results and to lighten the notation. Anytime we consider scattering we deal with plane 
waves, whether as states or as boundary conditions, which makes it more convenient 
to work with the momentum k instead of the energy E. 

Some of the poles of the ^-matrix can be identified with the bound states of the 
system: their locations in the fc-plane correspond to the bound state energies. The 
remaining poles turn out to be physically significant, as well. They correspond to the 
so-called anti-bound and resonance states, and hal^of the latter are sometimes called 
anti-resonance states. The poles will be the central objects of interest of this paper. 
We therefore provide a short summary of relevant results. 

Let us consider scattering by a potential that vanishes outside of a finite intervaQ 
[a, b], where a < 0, 6 > 0. The asymptotic wave function then has the form 

iPix) A±e*^'^ + 5^e-^^^ X ±oo, x i [a, 6]. (2) 

The map that relates the incoming (— ) and outgoing (+) coefficients is the S'-matrix: 

Now, let us define the Jost functions as the two independent solutions that behave 
asymptotically as e^*'^^', respectively, at ±oo. The nonzero Wronskian shows that 
0-1- (x, fc)* = 0±(x, — fc) is independent of 0±(a;, fc) and therefore we can express 0+ via 
0_ and its complex conjugate, and vice versa: 

<\)+{x,k) = a{k)(f)_{x,-k) I3{k)(j)^{x,k). (4) 

One can write a similar equation to express (p^ with different coefficients. Evaluating 
the Wronskians of each side with the appropriate choice of (j)±{x,±k), reveals that 
those coefficients can be expressed using a{k) and P{k) [3j. 

The components of S are related to the well known transmission and reflection 
coefficients (amplitudes) that are usually defined via the wave functions for a wave 
incident from the left or right: 

I ( \ — j '2*'''^ + -R+e"*''^ for X < a , , / x _ j T_e~*''^ for x < a , , . 

- < rj^^^ikx for X > 6 ; ^^^^^ ~ 1 e^^^^^' + i?_e*'=^' for x > 6 . 

Let us then consider a wave incident from the left. Notice that ipiix, k) = 0-(x, —k) + 
R+{k)(j)-{x, k) and ^/'^(x, fc) = T+{k)(j)+{x, k) satisfy the Schrodinger equation. More- 
over, expressing in terms of 0_ using ^ shows 

T+ = l/a{k), R+{k) = I3{k)/a{k). (6) 

Using as the wave function and proceeding just as we did with ipi reveals that the 
transmission coefficient T_ coincides with T+, so we define T{k) := T^{k) = T^{k). 
Also, in agreement with the conservation of probability, i?_ only differs from i?+ by 



^Resonances occur in pairs, at momenta ±fci — ik2, where fci, fc2 > 0. To distinguish the two poles 
of such a pair, those in the fourth (third) quadrant are known as (anti-)resonances. 
'^Then it clearly decreases 'fast enough' at infinity. 
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a phase and can be found, performing a similar calculation for a wave incident from 
the right. Using ip^ and ipR as the wave function in ^ and writing the corresponding 
values of the coefficients A± and B±, via ([s]), for each case allows us to express the 
^-matrix in the form 

'-{It)' 

a more familiar form, convenient for practical calculations. 

One can show (e.g. [20j) that the Jost functions are analytic functions of k in the 
upper-half fc-plane. The functions a{k) and P{k), therefore, are also holomorphic for 
Im(A;) > 0. It is clear now that the poles of the S-matrix are fully determined by the 
zeros of a{k). 

Now let us assume that a;(A;o)=0, for ko somewhere in the upper plane. From Q 
it follows that for k = k^ the Jost functions are linearly dependent. Now, recall that 
the Jost functions are solutions of the Schrodinger equation and since they are linearly 
dependent for k^ we have a solution that vanishes asymptotically at both infinities, 
i.e. a normalizable, bound state exists with energy k^. However, the Hamiltonian H 
is assumed to be self-adjoint, and therefore has only real eigenvalues. Since = 
1 + |/3(A;)p, a{k) ^ for any a; G M and thus the energy can only be real for ko purely 
imaginary. 

In addition a*{k) = a{—k) (from </)+(x, k)* = (f)+{x, —k)) and the Schwarz refiection 
principl^ imply that the analytic continuation of a{k) to the lower part of the complex 
plane is symmetric across the Im(A;)-axis. 

The complexification of the momentum, and therefore the energy, is more than a 
mathematical "trick", however. It has sound physical meaning. In the introduction we 
mentioned wave functions corresponding to complex eigenvalues and the corresponding 
Gamow vectors. The real parts of the complex eigenvalues z = E — iV /2 correspond to 
the physical energies, and their imaginary parts correspond to the decay widths, so that 
1/r is the lifetime of the decaying state [22j. The Gamow vectors satisfy Schrodinger's 
equation with outgoing boundary conditions (see, e.g. |23) and |22])- 

Our analysis so far does not imply anything about the existence of poles in the 
lower complex fc-plane. The only requirement is that the poles have to be symmetric 
with respect to refiection across the Im(fc)-axis. In particular, they can still lie on it, 
much like the poles associated with the bound states. If the potential is such that the 
S-matrix has poles in the lower plane, is there a pole-state correspondence as in the 
bound state case? 

As already mentioned, the states signalled by resonances are the Gamow vectors. 
This is easy to see in the simplest case of a square well. As Zavin and Moiseyev 
|18| show, outgoing boundary conditions lead to the same equations that determine 
the poles of the ^-matrix. In addition, the anti-bound states can be obtained using 
incoming wave boundary conditions. 

Let us point out that all poles except those on the positive Im(fc)-axis lead to non- 
normalizable wave functions. A way of dealing with the spatial divergence is considering 



®The Schwarz reflection principle states that for a holomorphic function in the upper complex 
plane, continuous and real valued on the real line, one can write an analytic continuation for the 
whole C plane such that f{z*) — f*{z). 
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the time dependent wave functions: 

^,(a;,t) = e-^^V^*/2^,(x), (8) 

which obviously decay with time, as long as F > 0. In terms of the wave number k 
we have F = — 4Re(A;)Im(/c)>0 only for Re(A;)>0, i.e. for resonance poles. The anti- 
resonances, with Re(fc)<0, correspond to the time-reversed behaviour of the resonances. 
The interpretation of the remaining divergent wave functions, those of the anti-bound 
states, with Re(fc)=0, remains obscure, though they are closely related to the bound 
states and resonances. In addition, as indicated by Nussenzveig even the resonance 
poles' interpretation breaks down for |F/i?| ^ 1. 

The complex eigenvalues do not contradict the well known reality of the eigenvalues 
of self-adjoint operators. The Hamiltonian, like any operator, is defined not just by 
its action but also by its domain. Strictly speaking then, two operators with the same 
action can have two different domains, which makes them two distinct operators. The 
domain for which we have determined the self-adjointness of the Hamiltonian, does 
not contain those complex energy states, since they are not normalizable in the usual 
sense. If we include the Gamow vectors in the domain the operator will no longer be 
self-adjoint. 

To conclude this preliminary section let us make a connection with another im- 
portant quantity - the resolvent. Resonances are also studied through its poles. The 
reason is that, while the S'-matrix does not exist for all Hamiltonians, the resolvent is 
a more general mathematical object and as such can be defined for much wider class 
of operators. And when scattering is possible, the poles of the S'-matrix coincide with 
those of the resolvent. 

Consider the time evolution of the state vector as determined by the propagator, the 
unitary operator of the form U{ti,t2) = e~*(*2-ti)-f^/'i for t-independent Hamiltonian 
H. Its Laplace transform defines the resolvent operator of the Hamiltonian, 

RH{E) = {H-Eir = J2\^^, (9) 

where / C M is a discrete or continuous interval depending if the energies are part of 
the discrete or continuous spectrum, respectively. Comparing ([T| and ^ shows that 
eigenvalues of the Hamiltonian generate poles of the resolvent. The resolvent, however, 
treats all the complex eigenvalues on an equal footing. 

The resolvent and the S'-matrix, as a function of A; = \/2mE/h, have the same 
singularities for all physical purposes, as we mentioned. This is to be expected since 
Rh{E) is related to the Green's function G{x, y): the Schwarz operator kernel of Rh{E) 
is the Green's function G{x,y), satisfying {H — E)G{x,y) = 6{x — y) (see, e.g. [24J). For 
scattering a natural assumption is that the initial state was the wave function defined 
in the far past, i.e. at t = — oo, and similarly the final state, at t = +oo. Thus the 
resolvent is also related to the scattering matrix via S = f/(— oo,+oo) and a Laplace 
transform. El 

For our purposes, however, finding the pole content of the scattering matrix is sim- 
pler than finding that of the resolvent. For piece- wise flat potentials it can be done by 

^For the mathematical rigour that we ignore, the reader may refer to the mathematical literature. 
Quite comprehensible lecture notes are available by Tang and Zworski [5S], for example. 
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simply matching the values of the wave function and its derivative at the discontinu- 
ities. Simple algebraic equations resuh, that can be handled easily by software. The 
resolvent makes the relationship between its poles and the spectrum clearer but the 
^'-matrix allows for easier calculations. 



3 Flows for the finite square well/ wall 

A particle moves freely along the real line except in [—a, a]; it interacts by the piece- wise 
flat potential 

n^) = l ~n^' ^^''''''Km. . (10) 
I 0, X G (-CX), -a) U (a, oo). ^ ' 

For positive/negative values of U we have a potential well/wall, as shown in Fig. [Tj 
The matching conditions for the wave function and its derivative are given by 

i){±a - 0) = i){±a + 0), i)\±a - 0) = ^'(±a + 0) . (11) 



We apply (11) to the wave function in the well t/^^ := Aexp(ii^x) + B exp{—iKx) and 
where K = \/k^~+~2U, which allows us to determine the wave function on the 
whole real line. This way we also determine the scattering matrix via ([T]) : 



Q= 1 2ikK 1 (]0) 

2kKcos{2Ka) - i{k^ + K^) sm{2Ka) y'^^'^J^f^^^"^ 1 >■ ^ ' 

The poles of the ^-matrix must obey the nonlinear algebraic equation 



kVk^ + 2Ucos{2aK) - i{k^ + U) siYi{2aK) = 0. (13) 

They will depend on a and U (and if we consider the asymmetric interval, also on b; 
see Sect. 4 below). In this section we consider the symmetric interval b = —a = 1.5 in 
dimensionless units. Fixing the width of the well, we will have a one-parameter family 



of solutions of (13). In the complex momentum plane it generates a flow, a collection 
of curves representing the positions of the poles as a function of the depth/height U, 
which we will discuss below. 

Using Mathematica, Maple or other similar software, it is easy to find all the sin- 
gularities kn- We used Mathematica's built-in functions to find all solutions in a 
bound region around the origin. Ideally, one can obtain the full flow by increasing 
the depth/height incrementally and generating a pole configuration at each step. The 
resulting points can then be plotted, displaying the flow. One can plot as many points 
as desired, making the plot more detailed or covering wider range of parameters. We 
introduce a depth/height cutoff in order for the algorithm to be finite. The cutoff is 
chosen so that the plots include the interesting features. At least in theory, one can 
decrease the step that changes the parameters to achieve very high level of detail. In 
practice, however, the average computer has the capability to generate only a limited 
number of configurations. For sufficient detail we plotted small portions of the flow at 
a time and then combined them. 
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Figure 3: Flow trajectories of the anti-resonance poles for the potential well/wall 
with variable depth/height. Only the trajectories for Iie{k) < are shown since there is a 
complete symmetry under reflection across the Im(k)-axis. The arrows indicate the direction 
of the flow for —U varying from —oo to oo (a deep well to a tall wall). 



The singularities of the S'-matrix are almost always simple pole^ with a typical 
configuration shown in Fig. [2j As expected from the general properties of the scat- 
tering matrix, discussed in the previous section, the singularities lie on the imaginary 
momentum clXlS clS well as in the lower complex plane. Following Nussenzveig [17J we 
plot the flow (i.e. the positions of the resonances in the complex plane) for depth/height 
U G (— oo, oo), treating the square potential well/wall, i.e. the well and wall together. 

First, let us consider the poles with nonzero real part, i.e. the resonances (or reso- 
nance poles/states). Recall that resonances occur in pairs, at two momenta with equal 
and negative imaginary parts, and equal and opposite real parts. Those in the third 
quadrant are known as anti-resonances. The fourth-quadrant poles are referred to sim- 
ply as resonances, when no confusion is likely, and sometimes as "physical resonances" 
when distinctions have to be clear. The flow is shown in Fig. [3] for the anti-resonances, 

^"For a discrete set of U values there is a double pole. 
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Figure 4: Flow trajectories of the bound and anti-bound states for a potential well 
(wall) with depth (height) -U>0 (U>0). The hollow points trace the positions of the 
branch points. The dashed line indicates the coalescence point k = —i/a. 
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Figure 5: Coalescence of S-matrix poles. Two anti-bound poles merge into a double pole, 
that consequently gives rise to a resonance + anti-resonance pair. 
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whose behaviour in time is the time-reversal of that of the resonances. These suffice 
for a complete description, because of the symmetry under reflection across the Im(A;)- 
axis. When the well is very deep the poles are located close to the horizontal line 
Im(A;) = — 1/a. As the well becomes shallower the poles move monotonously down- 
wards along the trajectories as indicated by the arrows in Fig. |3| In the case of a wall 
of increasing height, the poles travel upwards, cross Im(A;) = —1/a and diverge to the 
left. 

The horizontal line Im(A;) = —1/a bounds the pole trajectories for the well. More- 
over and as the depth increases a resonance + anti-resonance pair of poles travels to 
k = —i/a from each side of the Im(/c)-axis and coalesces to form a double pole, conse- 
quently creating an anti-bound pair, as in Fig. |5| This only happens for a countable 
set of values Un^n G Z, that can be found numerically. The point A; = —i/a is a coa- 
lescence point for the flow and it represents a transition of a pair of resonance states 
into a pair of anti-bound states. For a very deep well the imaginary part gets closer to 
the boundary, accounting for the characteristic flattening of the pole distribution for 
f/<0. 

The discrete subset of values of the depth/height f/„, for which a double pole exists, 
corresponds to each of the disjoint sections of the flow: Ui corresponds to the trajectory 
closest to the Im(A;)-axis (and its 4th quadrant counterpart), U2 to the second closest, 
etc. This is easy to see if we notice that once the resonance reaches the boundary 
and disappears there is no other resonance on that trajectory, as the whole resonance 
configuration just travels up and to the right along fixed branches of the fiow. 

Another way of labeling the trajectories is by their asymptotic behaviour. We find 



from (13) that the pole locations kn — ?■ ±7m/2a — ioo, n G N, as f/ — )■ 0. The rate 
at which the poles move also increases as f/ — )■ 0. In fact, as pointed out by Regge 
[21| . very small changes in the depth result in very large shifts in the locations of the 
resonances. Despite their sensitivity, however, they follow a smooth trajectory and 
while the fiow rate of the poles is large for shallower wells, the trajectory is well defined 
for all values of U. As — )■ 00, we verify that resonances are absent in the free 
particle case. 

For large n, the locations of the resonance + anti-resonance poles can be approxi- 
mated [3J . For a finite range constant potential we have Re A;„ ~ n and Im kn ~ —In n 
as n — )■ 00. It can also be shown that the distance between the asymptotic values of 
the branches of the fiow lim[/^oo(^n — ^n-i) is independent of n. 

As discussed, the parameter U can be taken into negative values to reproduce the 
square wall as a deformation of the square well. The fiow generated for f/ < comes 
from —ioo and shares vertical asymptotes with the fiow for f/ > 0. For that reason 
all the branches are parametrized by their asymptotic values, as before. The shape of 
the trajectories is shown in Fig. [s] In the infinite wall limit {U — )■ —00) we find the 
asymptotic values |A;„| — )■ 00. Clearly, for an "impenetrable" barrier, resonances do not 
exist. 

Since the fiow does not converge to a single point as in the well case, no bound 
and anti-bound states will be created by the resonances. For the attractive potential 
the resonance states are attracted towards each other to produce bound states. The 
location of the attractor (or coalescence point k = —i/a) is independent of the potential 
strength (at least for the cases we consider) but depends on its width. Conversely, the 
repulsive potential pushes the resonances away, preventing them from coalescing, and 
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so from subsequently producing a bound state. 

Now let us turn our attention to the singularities with Re(fc) = 0. They correspond 
to the bound states (Im(fc) > 0) and the anti-bound states (Im(/c) < 0). For a potential 
well, consider a typical bound-state-not the ground state, nor the first excited state. 
The corresponding pole location has (Im(A;) > 0), and as the well-depth decreases, 
(Im(/c) decreases monotonically. When (Im(A;) becomes negative, the state changes 
from bound to anti-bound. For a small range of U thereafter, it has the strange property 
of an "energy" that becomes more negative as the well becomes shallower. Then it 
meets another anti-bound state ascending from —ioc at the k = —i/a coalescence 
point. There a pair of anti-bound states "annihilate" giving birth to a resonance + 
anti-resonance pair of poles (Fig. |5]). The decrease in the number of bound states is 
tied to the decrease in the number of anti-bound states. The exceptions are the last 
two bound-state poles-one having as a limit and the other rapidly traveling to —zoo. 
For very shallow wells, a surviving bound state always exists. The last bound state 
becomes an anti-bound when the well turns into a wall, and heads to —ioo as the wall 
grows, with the Im(A;) axis as an asymptote. 

One can see that as the well becomes shallower the highest bound state will turn 
into a zero-energy state. Zero-energy states have infinite spread, i.e. their de Broglie 
wave length becomes infinite. They represent a critical point: the normalizable wave 



functions become divergent anti-bound states just above the well. From (13), we see 
that zero-energy states are admitted when U takes the values n^vr^/Sa^, n G N. 

Besides the bound and anti-bound state flows. Fig. |4]also depicts the trajectories of 
the branch points, distinguished by the hollow points. The branch cut is determined by 
i.e. Im(A;) e (-oo, -y/2U) U {V2U, oo). As we can see all the bound states 
are contained within the parabola, their energy cannot be smaller than the depth of the 
well, since the wave functions for E < Umin are not normalizable. As the anti-bound 
states are non-normalizable anyway, they are not affected by the branch cut, as the 
anti-bound states for — 1 < < illustrate in Fig. |4j An interesting observation is 
that the bound/anti-bound pairs are always contained within the branch points for all 
values of U . 

For the finite barrier, i.e. f/ < 0, only resonance (and anti-resonance) poles occur. 
Recall that when the resonance flow converges, a pair of trajectories collide from each 
side of the imaginary axis which results in the creation of a bound/anti-bound pair. 
For the potential wall the resonance flow diverges, so there are no bound states, and 
therefore, no anti-bound states either. 



4 Flows for the two-piece well+wall potential 

Now let us investigate how the presence of both attractive and repulsive parts affects 
the flow of the poles of the S'-matrix. Again, we consider the simplest case possible, a 
potential that has two constant non-zero pieces, the well+wall shown in Fig. [1} 

—U, X G [a, 0], 

V{x) = { H, xE (0,6], (14) 
0, X E (— oo, a) U (6, oo). 
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When either f/ or if change sign, the corresponding piece of the potential changes 
character, from a well to a wall, or vice-versa. Therefore, to illustrate the changes 
to the pole distribution, we will provide 2 graphs: Fig. 6 depicts the flow for the 
well+ well/wall, and Fig. 7 the wall+well/wall. 

The matching at the discontinuities provides the analytical expression for the S- 
matrix. Using the ansatz ([s]) we find the scattering matrix in the form of ([T]): 

= -e^'""^ [k{K'^ - k^) sm{aK) cos(6/€)+ (15) 
+ ik{K^ - K^) sm{aK) sm{bK) + K{K'^ - k^) cos{aK) sm{bK)] /T 
R = -ie-^'^^ [/t(A;2 - K^) sin(aK) cos(6/t)- 

- ik{K^ - K^) sm{aK) sm{bK) + K{k'^ - fc^) cos{aK) sin(6/«)]/r 
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Figure 7: Flow trajectories of the anti-resonance states. The flow is for a potential 
consisting of a fixed wall and variable well changing into a wall (fixed H , varying U ). 



with K = - 2if , K = + 2U and 

r = sm{bK) + k^)K cos{aK) + K{k^ + K^) sm{aK)] + (16) 

K cos{bK) [2kK cos{aK) + ^(A;^ + K^) sm{aK)] . 

The poles are then determined by the equation T = 0. The resonance flows are 
shown in Figs. [6] and [7j In the first case we have a potential that consists of a fixed 
well of depth U combined with a variable-depth {H < 0) well transitioning into a wall 
(H > 0). The second figure (Fig. [T]) shows the "reverse" case, where we fix the barrier 
and vary the well with depth —U < into a barrier with height —U > 0. 

Let us concentrate on the first case for now. The flows are now separated into 
two distinct families. For a fixed well + varied well we have the familiar behavior for 
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Figure 8: Resonance poles and the reflection coefficient. The well+wall potential 
has a well with fixed depth with the barrier increasing in height. Increasing the wall height 
"unzips" the line of resonances into two lines. 
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Figure 9: Resonance poles and the reflection coefficient. The well+wall potential has 
a well with fixed depth and a wall with fixed height, hut increasing width. The "unzipping" also 
occurs in this case. 
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Figure 10: Flow of the bound and anti-bound poles for the well+wall potential. The 

potential consists of a well with fixed depth and a well transforming into a wall with variable 
depth/height (fixed well + variable well/wall). 




Figure 11: Flow of the bound and anti-bound poles for the well+wall potential. The 

potential consists of a wall with fixed height and a well transforming into a wall with variable 
depth/height (fixed wall + variable well/wall). 
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Figure 12: Magnified flow for the 2-piece potential. Trajectories of the bound, and 
anti-bound states for the fixed well part (dashed line) and fixed wall (solid line) potentials 
around the origin. 



very deep (symmetric) well-all the poles tend to move towards k = —i/a, with the 
characteristic annihilation of resonances occurring as before. However the presence of 
the fixed well destroys the monotonous behavior of the flow and generates the oscillating 
pattern seen in both Fig. [6] and Fig. [7j The oscillating pattern shows that despite the 
fact that \U\ -C \H\, there is an observable influence of the finite part even when the 
variable part of the potential is of much greater magnitude. 

H ^ shows separation of the resonances into two families: rapidly diverging, due 
to the increasing wall, and slowly moving (finite well) ones. As we will confirm later, 
each family can be associated with one of the two pieces of the potential. The slowly 
moving half of the resonances need to recover the configuration for a finite well (of 1/2 
the width) so they remain finite for all H. The transition from H < to H > result 
in the local minima of Im(A;) observed in Fig. |6} 

As the well turns into a wall we observe rapid growth of new branches from —zoo. 
Those consist of poles that can be associated with the varied part of the potential. 
Unlike the symmetric well, however, the flow does not share asymptote with its < 
counterpart p] 

Again we observe two types of resonance poles. On one hand we observe that 
some of the poles follow diverging trajectories heading to the left with Re(A;) — t- — oo 
and lm{k) — )■ 0. The locations of the rest of the poles move upwards until they 
reach maxima and then descend with decreasing rate to reach a configuration that 
changes very slightly with the varied height of the barrier. However we observed 
similar behavior for the 1-piece potential barrier. This seems to suggest that indeed 

^"'^We have to truncate the vakies of H, in order for the algorithm to terminate, but we did not 
observe convergence between the H < and H > branches of the flow. 
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the poles lim.H->oo \k\ oo correspond to the barrier and the rest - to the welL Further 
confirmation can be found in Fig. |8]- it shows "snap shots" of the locations of the poles. 
The pole configuration exhibits a zipper-like behavior due to the separation of the poles 
along the two distinct flow branches, as discussed. The "top" poles are pushed away 
by the branch cut of \/k^ — 2H (introduced by the varied well/barrier) and does not 
affect the "bottom" poles. In addition Fig. |9] shows that widening the barrier, for fixed 
H does not affect the "bottom" poles while significantly increases the count of the "top" 
poles. 

The pole separation is quite surprising since the behaviour of the resonances is de- 
termined by a non-linear relation and yet they behave in a remarkably linear fashion. 
Only when the two pieces are of comparable size can the resonances not be distin- 
guished without information about their trajectories. However, labeling the poles will 
be incorrect since the poles corresponding to the fixed well and those corresponding 
to the varied well barrier can exchange roles and only make sense when \U\ ^ \H\. 
Indeed the poles that survive the limit if — t- are on the same trajectory as the poles 
that vanish in the limit — t- oo. 

Note that Fig. |9] indicates, as mentioned earlier, that the "top" resonances corre- 
spond to local minima of the reflection probability. Thus, the name resonance poles 
is justified-they are fully transmitted modes responsible for a transmission resonance. 
For low potential barriers the peaks begin to widen and overlap, as in the case of the 
"bottom" poles, gradually destroying the visual correspondence observed in Fig. [9j The 
last two figures also illustrate the physical significance of the resonance poles, as well 
as the branch cuts. The resonances indicate maxima in the transmission coefficients, 
while the brunch cuts act as an impenetrable barrier for the resonance states associated 
with the varied portion of the potential. 

Now let us consider the bound states. Fig. [TO] shows that when two wells interact 
the matching distorts the trajectories and they have inflexion points for certain values 



of the potential. This "complex plane" scattering is not present in the case of Fig. 11 
since the barrier does not have bound states. While there are also zero-energy states 
for a discrete set of potentials, they cannot be found in closed terms, due to the more 
complicated matching conditions. 

There is a noticeable difference around U = 0-instead of the least-bound state 
turning changing into an anti-bound one and then disappearing to infinity, we have 
an ascending anti-bound state annihilating it for a very small height of the wall. The 
most-bound state survives the transition from a well to a wall due to the presence of 
the fixed well to the left. 

The features around the zero potentials are interesting enough to be plotted in 
greater detail in Fig. [T2j Apparently, the absence of bound states associated with the 
barrier part leaves the familiar monotonous behavior of the bound/anti-bound states 
unchanged. 



5 Conclusion 

We studied the flow of ^-matrix poles in the plane of complex momentum for the 
simple potential(s) depicted in Fig. [l] Our results are summarized by the remaining 
Figures. 
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The most elementary potential considered was the square well/wall, and our Figs. 
[2]|5] recover the seminal results of Nussenzveig [T7]. Some of Nussenzveig's conclusions 
were also confirmed by less direct, graphical methods in jT8] (see also |13j). 

New is a similar study of the "two-piece" potential, consisting of two adjacent 
well/ walls, with strengths controlled by independent parameters (U and H in Fig. 
[T]). The flow of ^-matrix poles was calculated as a function of various potential pa- 
rameters, one at a time, and plotted in Figs. [6 12 



The Appendix also verifies that the "one-piece" and "two-piece" flows are consistent 
with the pole structure of a 5-function potential and a ^'-function potential, respec- 
tively, if the parameters are controlled in the appropriate way. The calculations are a 



simple, yet interesting check of our methods. Figs. p]3p5] depict the relevant flows. 

Unsurprisingly, the flows for the more complicated potentials are more complicated. 
Certain features can be understood simply, however, and we hope that our extension of 
the old analysis of Nussenzveig [17] will help lead to further insight that can be applied 
generally. 



Appendix: 6- and (^'-sequence potentials 

We can perform a simple, interesting check of our results by recovering the well- 
known textbook examples of the Dirac delta function (or 6) potential and its derivative, 
6', as limits of the "one-piece" and "two-piece" potentials, respectively. 

Consider first V{x) = —X6{x). We realize this 5-potential as the a — )■ limit of 
a (5-sequence of the "one-piece" potential: —U [6{x + a) — 9{x — a)], with the strength 
A = 2aU kept fixed. To see what we should find, we solve the Schrodinger equation for 
X G IR/{0} and demand that the wave function is square-integrable to show if) oc e""**'^'. 
Integrating equation ([T| over an infinitesimal interval (—£,£) around yields 

V^'(e) - V^'(-e) = -2mAV^(0), (17) 

which determines the energy E = — of the bound state. A straightforward 
calculation shows the transmission and reflection amplitudes 

t = (1 - m\/ik)-\ r = {ik/m\ - 1)"^ (18) 

have a pole only at the bound state energy, i.e. no resonances exist for the (5-potential. 
This is also confirmed by the fact that the probabilities for transmission |tp and re- 
flection |rp are monotonous functions, i.e. there are no transmission maxima. 

The graphical representation of the flow leads to the same conclusions. The pole 
spectrum of the square well flows asymptotically to that of the 5-potential. Setting 
all = 1 and taking the limit f/ — )■ oo produces the trajectories of the resonance poles 



shown on the left side of Fig. 15 The resonances approach complex infinity as the 



depth increases and no finite resonances exist in the 5-limit. Similarly, a single bound 



state survives (Fig. 13) in agreement with the above discussion. The final anti-bound 



state rapidly diverges to —ioo leaving no trace. 
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Figure 13: Flow of the bound/ anti-bound states of the 5 -sequence potential. A 

single hound state survives as the potential approaches the delta-function limit. 



Im(k) 




Figure 14: Flow of the bound state of the 6' -sequence potential. A single bound state 
survives; however, it is a zero- energy state. 
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Now let us consider the derivative of the 5-potential as arising from 

- U [e{x + a) - 2e{x) + e{x - a)] , (19) 

a special case of the "two-piece potential". We can show that the Schrodinger equation 
with (^'-potential has the same solutions as for the (5-potential but the matching has 
to be done differently. Assuming the wave function is continuous at x = 0, we once 
again integrate the Schrodinger equation around zero in order to obtain a matching 
condition for the derivative: 

^'(-e) = V^'(e). (20) 

This is quite different from the 5-potential, since the matching condition demands 
a wave function that is also smooth at the origin, which is only possible for k = 
0. Matching for a plane wave incident from the left and right shows the absence of 
resonances and anti-resonances as well. 

While somewhat trivial, this potential illustrates another zero-feature limit, that 



of the well+wall. The flows generated by the limit are shown in Fig. 14 for the 



bound/anti-bound poles and Fig. 15 (together with the 6 case), in agreement with the 
analytic arguments. 
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